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AN  EXACT  ANALYTICAL  EXPRESSION  OF  THE  SHEAR  COEFFICIENT 
IN  THE  MINDLIN  PLATE  EQUATION 

INTRODUCTION 


The  Mindlin  plate  equation  is  a  model  of  a  thin  plate  with  the  addition  of  rotary  inertia  and 
shear  effects  on  the  dynamics  of  the  plate  response.1  Incorporated  in  the  Mindlin  plate  equation 
is  a  shear  coefficient,  which  is  an  adjustment  parameter  in  thick  plate,  beam,  and  shell  equations 
of  motion  that  is  included  to  compensate  for  stress  distribution  in  the  cross-sectional  shape  of  the 
object.  It  is  sometimes  designated  the  Timoshenko2,3  shear  coefficient,  a  concept  introduced  by 
Timoshenko  to  account  for  the  non-uniform  shear  stress  distribution  in  a  beam.  Timoshenko 
estimated  the  shear  coefficient  as  a  function  of  Poisson’s  ratio  in  the  beam.3  Cowper4  also 
estimated  the  shear  coefficient  as  a  function  of  Poisson’s  ratio,  although  he  derived  a  slightly 
different  expression.  Over  the  years,  it  has  become  evident  in  beam  theory  that  the  shear 
coefficient  is  theoretically  dependent  on  more  than  Poisson’s  ratio.5'8  The  problem  of 
determining  the  shear  coefficient  in  a  Mindlin  plate  equation  has  been  addressed  by  Stephen,9 
who  determined  what  he  called  a  “best  shear  coefficient”  based  on  matching  a  mode  of  the 
Mindlin  plate  theory  to  the  exact  Rayleigh-Lamb  frequency  equation  for  the  flexural  wave 
response. 

The  above  theories  are  all  based  on  excitation  of  the  flexural  wave  in  a  structure,  and  the 
corresponding  shear  coefficients  are  determined  at  the  specific  wavenumber  and  frequency  of  the 
flexural  wave.  They  do  not  account  for  the  shear  coefficient  as  a  function  of  all  wavenumbers 
and  frequencies,  which  is  a  response  and  excitation  condition  that  exists  in  structures  that  are 
loaded  by  turbulent  boundary  layers  or  acoustical  forces. 

This  report  derives  an  exact  analytical  expression  of  the  shear  coefficient  in  the  Mindlin 
plate  equation  subjected  to  planewave  excitation  at  any  wavenumber  and  frequency.  This  is 
accomplished  by  computing  the  displacement  field  of  the  plate  using  the  Mindlin  plate  equation 
and  the  thick  plate  equation,  and  then  setting  them  equal  to  each  other.  Because  the  shear 
coefficient  is  explicit  in  the  Mindlin  plate  equation  and  implicit  in  the  thick  plate  equation,  it  can 
be  solved  for  as  a  function  of  wavenumber,  frequency,  and  plate  parameters.  A  numerical 
example  is  included  to  depict  the  dependence  of  the  shear  coefficient  on  wavenumber  and 
frequency.  It  is  shown  that  the  shear  coefficient  is  only  slightly  dependent  on  frequency  and 
extremely  dependent  on  wavenumber  of  excitation.  Comparisons  of  previous  analytical 
expressions  are  also  included  in  the  numerical  example  to  illustrate  how  other  theories  compare 
to  the  one  derived  here. 
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SYSTEM  MODELS  AND  THEORETICAL  SHEAR  COEFFICIENT 


Two  system  models  are  developed;  one  contains  the  shear  coefficient  explicitly  and  the 
other  contains  the  shear  coefficient  implicitly.  The  first  system  model  is  that  of  a  Mindlin  plate 
whose  governing  equation  is10 


DV2  - 


p/?3  d2 
12  dt 2 


\u(x,t)  +  ph 


d2u(x,t) 
dt 2 


'  DV2  | 

kph  Ylkp  dt1 


\f(x,t),  (1) 


where  k  is  the  shear  coefficient,  u(x,t)  is  the  displacement  of  the  plate  in  the  z-direction,/x,0  is 
the  force  distribution  on  the  plate,  p  is  the  density,  p  is  the  shear  modulus,  h  is  the  thickness,  x  is 
spatial  location,  t  is  time,  V  is  the  spatial  gradient  operator,  and  D  is  equal  to 


Eh 


12(1  -v2)’ 


(2) 


where  E  is  Young’s  modulus,  and  v  is  Poisson’s  ratio.  The  system  is  modeled  as  infinitely  long 
with  a  continuous  forcing  function  varying  in  time  and  space;  thus,  the  displacement  and  forcing 
function  terms  are  written  as 

w(x,  0  =  U(p,  co )  exp  (icot)  exp  (ipx) ,  (3. 


and 

f(x,t)  =  F(p,  co)  exp(/  o)t)exp(ipx) , 

where  co  is  angular  frequency  and  p  is  wavenumber  with  respect  to  the  x-axis.  Solving  the 
transfer  function  response  of  displacement  divided  by  input  force  yields 

Ujppo) _ -12/M-12  Dp1  +  ? 

F(p,co)  (12 Dphp4  -pph*co2p2  -\2pph2co2)k  -UDphco2 p2  +p2h4co 4 


U_ (p, co)  _  ak+]> 
F(p,co)  ck  +  d ’ 

where 

a  =  -12  ph. 


(6) 


(7) 
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(8) 


b  =  -\2Dp2  +  ph3o)2, 


c  =  12  Dphp  4  -  pph4 a>2 p2  -12  pph2 co2  , 


(9) 


and 


J  =  -12 Dphco2p2  +  p2h4co 4  . 


(10) 


The  second  model  is  derived  from  the  equations  of  motion11  of  a  solid  medium,  governed  by 


<?2u 


//V  u  +  (A  +  p)VV  •  u  =  pz-r, 

d\L 


(11) 


where  X  and  p  are  the  Lame  constants,  •  denotes  a  vector  dot  product,  and  u  is  the  Cartesian 
coordinate  displacement  vector  of  the  plate.  Assuming  harmonic  response  in  space  and  time, 
equation  (11)  can  be  manipulated;  the  resulting  expression12  is  the  displacement  in  the  ^-direction 
at  the  middle  location  of  the  plate  divided  by  the  incident  force  and  is  written  as 

^Z(P> -^/2>^)  Q(/7, -h/2,o>)  /j2) 

F{p,(o)  pk(p,co) 


where 

Q>(kx ,  - h 1 2,a>)  =  %a2 ft p2 (P2  -  p2) sin(ah / 2) cos2 (Ph / 2) 

+  4 ap2(p2  -  p2)2  sin(/?/j /  2) cos2  (ah /  2) 

+  1 6a2  Pp4  sin(ah / 2)cos(a/z / 2)cos(/?/i / 2) 

+  2a(/?2  -  p2)3  sin(ph/2)cos(ah/2)cos(ph/2),  (13) 


and 


A(/?, 6>)  =  -  8a/? p2(P2  ~  p2)2 [cosfah) cos(/? h)  - 1] 

+  [(/?2  - p2)4  +  \6a2 p2 p4 ]sin(ah)sin(Ph).  ^4^ 

In  equations  (13)  and  (14),  cr  is  the  modified  wavenumber  associated  with  the  dilatational 
wave  and  is  expressed  as 

»  (15) 
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where  kd is  the  dilatational  wavenumber  equal  to  col cd,  where  cd  is  the  dilatational  wavespeed 
(m/sec);  is  the  modified  wavenumber  (rad/m)  associated  with  the  shear  wave  and  is  expressed  as 

H^-p2  •  <16) 

where  ks  is  the  shear  wavenumber  (rad/m)  equal  to  a>/cs,  where  cs  is  the  shear  wavespeed 
(m/sec).  The  relationship  between  the  wavespeeds  ( cd  and  cs)  and  the  plate’s  Lame  constants 
(A.  and  /u)  is  determined  by 


It  is  noted  that  this  shear  coefficient  is  not  only  a  function  of  wavenumber  and  frequency, 
but  also  of  Young’s  modulus,  shear  modulus,  Poisson’s  ratio,  and  the  thickness  and  density  of 
the  plate. 
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NUMERICAL  EXAMPLE 


A  numerical  example  is  now  analyzed  to  investigate  the  behavior  of  the  shear  coefficient 
calculated  using  equations  (19)  and  (20).  The  parameters  of  the  plate  are  listed  in  table  1.  Figure  1 
is  a  plot  of  the  shear  coefficient  versus  wavenumber  and  frequency  as  determined  using  equation 
(19).  The  parabolic  line  on  the  plot  is  the  flexural  wavenumber  location  calculated  by  finding  the 
maximum  value  of  the  displacement  in  wavenumber  at  each  analysis  frequency.  The  weak  line 
originating  at  the  origin  and  ending  at  f=  10,000  Hz  and p  =  75.4  rad/m  is  the  plate  wave 
dynamics  propagating  through  the  analysis.  This  plate  wave  wavenumber  can  be  predicted  by 


Pp=<°  1 


pa- 


-v2) 


(21) 


The  thick  plate  equation  of  motion  contains  the  plate  wave  dynamics  while  the  Mindlin  plate 
equation  does  not;  thus,  accounting  for  a  modeling  difference  around  the  plate  wavenumber.  The 
result  of  this  model  mismatch  is  that  the  theoretical  shear  coefficient  factor  will  not  be  accurate  in 
the  region  around  the  plate  wavenumber.  Figure  2  is  a  cut  of  figure  1  in  wavenumber  at  4000  Hz, 
which  illustrates  the  dependence  of  the  shear  coefficient  on  wavenumber.  The  square  marker  (□) 
is  the  plate  wavenumber  and  the  round  marker  (O)  is  the  flexural  wavenumber.  Figure  3  is  a  cut 
of  figure  1  in  frequency  at  the  flexural  wavenumber.  The  shear  coefficient  estimate  derived  by 
Timoshenko3  is  indicated  by  the  symbol  x  and  is  given  by  the  wavenumber  and  frequency- 
independent  equation 


,  _  5  +  5 v 
_  6  +  5t> 

The  star  symbol  (*)  is  the  shear  coefficient  derived  by  Stephen9  and  is  equal  to 


(22) 


k  = 


5 

6 -o’ 


(23) 


and  the  plus  symbol  (+)  is  the  shear  coefficient  estimate  derived  by  Cowper4  and  is  written  as 


10  +  10  o 
12  +  11  v ' 


(24) 


These  shear  coefficients  are  shown  at  low  frequency  because  they  were  derived  using  long 
wavelength  assumptions.  It  is  clear  that  these  estimates  are  close  to  the  exact  solution  derived  in 
equation  (20)  for  wave  propagation  at  (or  near)  the  flexural  wave  propagation  wavenumber  and 
frequency.  They  do  not  account  for  wave  propagation  at  other  wavenumbers. 
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Wavenumber  (rad/m) 


Table  1.  Plate  Parameters  Used  for  Numerical  Example 


Parameter 

Value 

Thickness,  h 

0.01  m 

Young’s  Modulus,  E 

7.0e8  N/m2 

Shear  Modulus,  G 

2.5e8  N/m2 

Poisson’s  Ratio,  v 

0.4  (dimensionless) 

Density,  p 

1200  kg/m3 

Frequency  (Hz) 

Figure  1.  Shear  Coefficient  Versus  Wavenumber  and  Frequency  for  the  Numerical  Example 
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Figure  2.  Shear  Coefficient  Versus  Wavenumber  at  4000  Hz  for  the  Numerical  Example 


Frequency  (Hz) 


Figure  3.  Shear  Coefficient  Versus  Frequency  at  Flexural  Wave  Resonance 

for  the  Numerical  Example 
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The  most  interesting  aspect  of  this  analysis  is  that  the  shear  coefficient  is  highly  dependent 
on  the  wavenumber.  This  dependency  was  previously  unknown;  however,  the  fact  that  the 
dynamics  of  the  system  are  dependent  on  the  wavenumber  of  the  excitation  is  not  a  new  concept. 
In  this  type  of  analysis,  low  wavenumber  excitation  generally  results  in  a  structural  response  that 
is  composed  of  primarily  dilatational  waves  that  contain  the  majority  of  the  energy.  (At  zero 
wavenumber,  the  response  is  entirely  a  dilatational  wave  response.)  Shear  effects  are  secondary 
until  the  wavenumber  of  the  excitation  becomes  moderate.  Previous  work  has  typically 
consisted  of  exciting  the  flexural  wave  in  a  beam,  measuring  response,  and  then  back-calculating 
the  shear  correction  factor.  It  has  not  consisted  of  exciting  the  structure  at  all  wavenumbers  and 
then  determining  the  correction  factor.  (It  should  be  noted  that  broadband  wavenumber 
excitation  of  a  structure  is  not  easy  to  implement.) 


CONCLUSIONS 


The  theoretical  shear  coefficient  for  a  Mindlin  plate  has  been  derived  as  an  analytical 
expression.  It  is  shown  that  this  term  is  dependent  on  wavenumber,  frequency,  Young’s 
modulus,  shear  modulus,  Poisson’s  ratio,  density,  and  thickness  of  the  plate.  Numerical 
simulations  showed  that  the  shear  coefficient  is  extremely  dependent  on  wavenumber  but  only 
slightly  dependent  on  frequency.  Previous  shear  coefficient  expressions  are  close  to  the 
analytical  expression  derived  here  around  the  flexural  wave  wavenumber  of  the  plate. 
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